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Abstract 

Explicit formulas for the rotation frequency and the long-wavenumber diffu- 
sion coefficients of global spirals with m arms in Rayleigh-Benard convection 
are obtained. Global spirals and parallel rolls share exactly the same Eckhaus, 
zigzag and skewed- varicose instability boundaries. Global spirals seem not to 
have a characteristic frequency uj m or a typical size R m , but their product 
u> m R m is a constant under given experimental conditions. The ratio Ri/Rj 
of the radii of any two dislocations (R4, Rj) inside a multi-armed spiral is 
also predicted to be constant. Some of these results have been tested by our 
numerical work. 

PACS numbers: 47.20.Bp, 47.54.+r, 47.20.Lz, 47.27.Te 
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Global spirals and spiral defect chaos (SDC) as intrinsic patterns have been experimen- 
tally observed recently in Rayleigh-Benard convection (RBC) |3],[|[T(]]. These observations 
were rather surprising because they were carried out in a parametric region where the fa- 
miliar parallel-roll pattern should be stable |J. An explanation for the unexpected presence 
of global spirals or SDC, in place of parallel rolls, is still to be provided. So far theoretical 
attempts on understanding these intriguing patterns rely heavily on numerical solutions of 
either the generalized Swift-Hohenberg (GSH) model HBH^,^| or the truncated Navier-Stokes 
equations governing the fluid dynamics ||. Although these numerical studies have repro- 
duced experimental results qualitatively and quantitatively, very limited theoretical insights 
have been obtained. While the formation of SDC in the system has since received con- 
siderable attention |8|||, little effort has been given to determining the essential properties 
of a single spiral. It is far from clear whether a global spiral has a characteristic rotation 
frequency u m or a typical size R m @. The knowledge of those properties, we believe, is 
necessary in order to describe the much more complicated SDC. 

Recently major progress was made by Cross and Tu (CT) in this front []12|| . Applying the 
phase dynamics method developed earlier in studying pattern formations in non-equilibrium 
systems |I3yT5|] , CT considered the dynamics of a spiral as the balance of two competitive 
motions: a radial phase-drifting of the rolls and an azimuthal climbing of the dislocation 
||12|| . CT's results also imply that uj m R m is a constant under given experimental conditions, 
but they did not give an explicit expression for u m R m . Furthermore, CT demonstrated 
qualitatively that the rotation frequency u m of a spiral is not directly related to mean flow, 
which is induced by distortions of the convective rolls [|13|], although mean flow is necessary 
for the formation of the rotating spiral 



In this paper we focus on dynamical properties of global spirals. We extend CT's results 
in two respects: We first make a one-mode approximation for spirals JIB] which leads to an 
explicit expression for uj m R m ; we also separate the phase fluctuations from the stable phase- 
drifting and calculate the phase diffusion coefficients. We test some of these formulas by 
our numerical solutions. Our results make it possible to discuss the dynamical properties of 



global spirals in detail. We predict that, inside a stable multi-armed global spiral, the ratio 
of the radii of two dislocations is a constant under given experimental conditions. We also 
find that global spirals, concentric rings and parallel rolls have exactly the same Eckhaus, 
zigzag and skewed-varicose instability boundaries. Presumably there is a competition among 
the various attractors corresponding to these states, the nature of which requires further 
theoretical study. 

To be concrete, we base our calculations on the two-dimensional generalized Swift- 
Hohenberg (GSH) model for RBC [13[], which has been proven very successful in charac- 
terizing convective patterns under quite broad conditions Numerical solutions of the 
GSH model not only reproduce both global spirals and SDC but also resemble experimen- 
tal results reasonably well In this model, the convective patterns are determined 
completely by an order parameter ip(r,t) in two-dimensional space r, which satisfies []13[ 



dtiff + U • W> = e - (V 2 + l) 2 V - # + £ 3 (W>) 2 V 2 ^ 
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where U is the mean flow velocity given by U = VC x e z while [I3| 
~d t - <x(V 2 - c 2 )] V 2 C = g m e z ■ [V(V 2 ^) x 
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In the GSH equations, the reduced Rayleigh number e = 2.7824e oxpt , where e exp t = 
(Ra/Ra c ) — 1 is the control parameter |6]f|, in which Ra and Ra c are the Rayleigh number 
and its critical value at onset [lj. Other parameters g, g$, c 2 , g m and the fluid Prandtl 
number a model the properties of the system and are all non-negative. For simplicity, we 
only consider Oberbeck-Boussinesq fluids here [2(J . 

A one-mode approximation has been used in studying spirals in chemical reaction- 
diffusion systems We now apply the same approximation for stable global spirals 
in RBC. Using polar coordinates (r, if), one can approximate a global-spiral solution by 
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where 8 m = k m (r)r + mip — \m\uo m t with m the number of spiral arms near the core. In 
general, the amplitude A m (r) and the wave number k m (r) should depend on r. The rotation 



frequency u m , however, must be independent of r for a stable spiral. Here we adopt the 
following conventions: k m (r) > 0, m > 0(< 0) if the spiral is right (left) -handed and uj m > 
0(< 0) if the spiral rotates in the same (opposite) direction of chirality. Experiments PJTO 



showed that a multi-armed spiral usually has dislocations at different radii. Across each 
of these radii, the number of arms decreases or increases (depending on whether m > 
or m < 0) by the number of dislocations on the corresponding radius. For mathematical 
simplicity, we denote each of these radii, for example R4, by the number of non-terminating 
arms, say i, in its inside vicinity: See Fig. |l| for a m = 3 global spiral. Then a multi-armed 
spiral with dislocations at different radii can also be decribed by the above solution, provided 
that m is replaced by % for each corresponding region Rj < r < B4. The amplitude and the 
wave number should be continuous across each boundary r = The frequency, on the 
hand, must be a constant for all regions. For r > Ri, only concentric rings (a target state) 
exist which corresponds to i = 0. With this understanding, our results below can be directly 
extended, by replacing m with every possible i, to multi- armed spirals found in experiments. 
Phase equation TT^JTH] describes slow variations of convective rolls from their perfect 



pattern. Assuming R m 3> 1, i.e., 77^ = 1/R m <C 1 for a global spiral, one can then introduce 
slow scales R = rj^r, $ = if and T = rj^t, and a slowphase variable 6 m (R, T) = ^0 TO (r, t) + 
\m\uj m R m T . Here the phase fluctuations m (R, T) have been explicitly separated from the 
stable phase-drifting \m\uj m R m T . A similar separation has been used in studying chemical 
waves with steady velocity Now the local wave number can be defined as |L4] 



q m = V r # m (r, t) = V R 6 m (R, T) = q m e r + 0{rj 2 m ). (4) 
Inserting the spiral solution (§) into the GSH equation (JJ), one may then match the result 



to order rj m . Since the calculations are carried out in exactly the same way as in Ref. |jT3 
we skip details here but only write down the final results. From the zeroth order of r/ m , one 
finds that the amplitude is slaved by the wave number and is given by [[[£ 



From the second order of T] m , one obtains essentially the phase equation 
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where 



B(q m ) = 2(1 -q 2 J--g 3 q 2 m \A m \ 2 , 



(7) 



while U' = VR,(C)e m x e 2 with (• • -)g m for the phase average in 8 m and, with g' m = g m /c 2 a, 



V R (C)e m = -7;9' m e z ■ q m x V R Vr ■ (q m |A 



(8) 



Although one may, in principle, convert Eq. (g) into the standard form as in Ref. [14], there 
is no need to do so here. 

To express the phase equation in the more familiar diffusion equation form, one needs to 
project the gradient operator V into local coordinates. In the light of Eq. (^), this can be 
easily achieved by d\\ = Or and d± = -R -1 <9$ ~ 9$ (since R ~ 1 for r ~ R m where the phase 
equation is valid). Again recalling Eq. (H), one finds immediately that Vr • q m = V R m . 
In Cartesian coordinates, this simply gives d 2 m + d\Q m which contributes to the diffusion 
of phase fluctuations. But in a polar coordinate system, an additional i? _1 9||0 m ~ q m term 
is also present which, however, contributes to the stable phase-drifting \m\uo m R m T\ From 
Eq. (H), one gets via this -R -1 <9||0 m ~ q m term that JL9 



\m\iO m R m = q m B(q m ), 



(9) 



whose corrections are of order 77^. Apparently this frequency of rotation is generated by the 
curvature of convective rolls although it appears independent of mean flow. For m = 0, it 



naturally leads to the "wave number selection" go — Qf with B(qj) = Jl^JT5 



Further algebra reduces Eq. m) to a phase diffusion equation [0,1/ 



d T Q m = D\\(q m )diQ m + D ± (q m )diQ m + D x (q m ) Vr^^O 



(10) 
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in which the diffusion coefficients are [19 
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(11) 
(12) 
(13) 



2 ^m dq 

This diffusion equation describes three types of long- wave- length fluctuations [0,0]: Eck- 
haus (Du), zigzag (D±) and skewed- varicose (D sv ) with, for D x < 0, 



(D ll -D ± ) 2 + 2(D ll +D ± )D x +D 2 x /AD X . 



(14) 



When all these diffusion coefficients are positive, global spirals are stable; but when any one 
of them becomes negative, spirals lose their stability against the corresponding fluctuations. 

One striking feature of this diffusion equation is that all the functions D», and D sv 
are independent of m. Indeed they all agree with those of parallel rolls. This means that 
global spirals, concentric rings and parallel rolls share exactly the same Eckhaus, zigzag 
and skewed-varicose instability boundaries. So the fact that the measured wave numbers of 
stable spirals are inside the stable region of parallel rolls 0] is not surprising but necessary. 
Considering that the system is non-potential and that the only difference between spirals of 
different number of arms is q m , a theoretical understanding of spiral-to-target or spiral-to- 
spiral transitions could be subtle. There are, however, some shortcomings in our analysis. 
The core instability [0] is omitted here. Short-wavelength fluctuations [p|,|l7| have also been 
neglected, but these seem irrelevant in the transition between global spirals of different value 
of m @,|10j- I n an y event, it is clear that further theoretical work on this issue is necessary. 

We now concentrate on the rotation frequency of the global spirals. One sees easily from 
Eq. (Q) that uj m R m is a constant under given experimental conditions. This has been implied 
by CT [|12j and verified by experiments [|TTJ and our numerical solutions: see Table |. The 
quantitative value of uj m R m apparently depends on q m . However, except for m — 0, one 
cannot determine q m merely from the phase equation. This problem has been emphasized 
by CT [ l2| . They found that it is essential to take the intrinsic defect of the spiral, i.e., the 
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dislocation defect, into account. A spiral can be stable only if the phase-drifting of the rolls 
is balanced by the climbing of the dislocation, i.e., ui m R m = v d (q m ), where v d (q m ) is the 
climbing velocity of the dislocation. This, together with Eq. (|9|), selects the wave number 
for a stable spiral. Unfortunately, an explicit analytic formula of Vd{q m ) seems intractable 
|TB| ]. Thus an accurate evaluation of q m seems beyond reach. 

One may still make the following approximations for small e [ |T^] , 



\m\u m R m ~ a(q f - q m ), v d (q m ) ~ f3{q m - q d ), (15) 

where a, [3 > and q d is defined by v d (q d ) = 0. For g 3 = 0, we get from Eq. @ qj = 1 and 
a = 4, while for g 3 > 0, we get qf ~ 1 — 7/e and a ~ 4(1 + ae) with 7^ = 3g 3 /4(3g + g 3 ) 
and a = 67/ + y7/. Similarly, one may write q d ~ 1 — 7^6 and /3 ~ /5o(l + Under 
these approximations, the selected wave number of a stable spiral is fully determined by 
UmRm = v d , which gives, with 7 m = (47/ + \m\(3^ d )/(A + |m|/3 ), 

q m ^ (aq f + \m\(3q d )/ (a + \m\(3) w 1 - 7 m e. (16) 

Then, from Eq. (|J), one finds for m^O that 

^*-rt i )K mv-^ _ (17) 

a + |m|p 4 + |m|p 

For a multi-armed spiral with dislocations at different radii, since the frequency of rotation is 
a constant over the whole spiral, the ratio Ri/Rj — (a+ \j\f3)/(a+ \i\0) is fixed under given 
experimental conditions. (The definition of Ri and Rj is given below Eq. (|3|).) This ratio 
depends on only two experimentally measurable quantities a and (3 ||11|| , which provides a 
strong test of our theory. 

To test the validity of Eq. (|]), we have solved Eqs. ([I]) and (Q) numerically by the same 
method described in Ref. 0, except that we set up a global spiral as the initial condition. 
We use g = 1.0, g 3 = 0, a = 1.0, c 2 = 2.0 and g m = 10 for our numerical solutions. 
We use mesh points iVxiV = 512x512 and grid size Ax = Ay = ir/8 for aspect ratio 
(radius/thickness) T = 32. For simplicity, we set all radii of dislocations inside a multi- 
armed spiral to be equal. For e = 0.3, a one-armed spiral shrinks rapidly into concentric 



rings. But for e = 0.5, one-armed, two-armed and three-armed spirals all stabilize with three 
given sizes R m /L = 0.45, 0.55, 0.65 where L = 32ir is the radius of the cell. We then measure 
the rotation frequency and, via Eq. (^), calculate the corresponding wave number q m . The 
results are listed in Table |. Although R m varies by about 50%, the product uj m R m is found 
to be a constant within 5%. Also from Eq. ( |TTD and the measured values of uj m R m , we find 
that (3 = 0.73 ± 0.08 and qd = 0.880 ± 0.001. In comparison, one gets qd ~ 0.905 from a 



formula in Ref. ||16|| . Furthermore, although global spirals are set up with three different 
sizes, they are all found to be stable within the time of our computer simulation, which runs 
about 2000 vertical diffusion times for each case. So it seems that global spirals do not have 
a typical size. Consequently, the frequency of rotation may also be non-unique. 

Finally we make a comment on the role played by mean flow in the dynamics of spirals. 



Evidently Eq. (||) does not explicitly depend on g m . This has been first observed by CT |T2| 



who hence assign a "secondary" role to mean flow. Nevertheless mean flow plays a subtle 
role in determining the selected wave number ( |16D for stable spirals. Indeed the value of 
qj — qd in Eq. (|T7]) is very sensitive to g m ]T5|] . So 0J m R m also has a sensitive g m dependence. 



Furthermore, one must have D± > for stable spirals. This, from Eqs. (|T^) and ([]), gives 
a constraint on the allowed frequency. For g m = 0, only u m R m < is permitted. Now if 
Lo m R m > is also necessary to avoid "unwinding" |12j| , stable spirals must be stationary. 



Recalling Eq. (|17D, a stationary spiral is possible only if q$ = q d . For #3 = and g m = 0, 
the relation qj = q d indeed holds [|TIJ and a stationary global spiral has been found ||. 



But for more realistic g% > 0, these two wave numbers are in general unequal. So a finite 
g' m is needed to observe any stable spiral, which might suggest why a low Prandtl number 
(g' m ~ 1/a) is necessary in experiments PJ5|JT0[|. 

In summary, we have calculated the rotation frequency and the long-wavelength diffusion 
coefficients of global spirals. We find that global spirals have exactly the same Eckhaus, 
zigzag and skewed-varicose instability boundaries as parallel rolls and concentric rings. So 
a transition between these patterns presumably involves a competition among their various 
attractors. Although global spirals seem not to have a characteristic frequency or a typical 



size, the product of them is a constant under given experimental conditions. The ratio of 
the radii of any two dislocations inside a multi-armed spiral is also predicted to be constant. 
Some of these results have been tested by our numerical solutions. Nevertheless, to fully 
understand the intriguing global-spiral pattern, an analysis of the core instability and a 
theory describing the spiral-to-target transition will be necessary. 

XJL and JDG are grateful to the National Science Foundation for support (under Grant 
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at the Pittsburgh Supercomputing Center and Cray-YMP8 at Ohio Supercomputer Center. 



9 



REFERENCES 

[1] M. C. Cross and P. C. Hohenberg, Rev. Mod. Phys. 65, 851 (1993). 

[2] G. Ahlers, Over Two Decades of Pattern Formation, A Personal Perspective [preprint]. 

[3] F. H. Busse, Rep. Prog. Phys. 41, 1929 (1978). 

[4] E. Bodenschatz et al, Phys. Rev. Lett. 67, 3078 (1991); E. Bodenschatz et al, Physica 
D 61, 77 (1992). 

[5] S. W. Morris et al, Phys. Rev. Lett. 71, 2026 (1993). 

[6] H.-W. Xi, J. Vihals and J. D. Gunton, Phys. Rev. A 46, R4483 (1992); H.-W. Xi, 
J. D. Gunton and J. Vihals, Phys. Rev. E 47, R2987 (1993). 

[7] M. Bestehorn et al, Phys. Lett. A 174, 48 (1993). 

[8] H.-W. Xi, J. D. Gunton and J. Vihals, Phys. Rev. Lett. 71, 2030 (1993); H.-W. Xi and 
J. D. Gunton, Phys. Rev. E 52, 4963 (1995). 

[9] W. Decker, W. Pesch and A. Weber, Phys. Rev. Lett. 73, 648 (1994). 

[10] B. B. Plapp and E. Bodenschatz [preprint]. 

[11] E. Bodenschatz and B. B. Plapp, Bull. Am. Phys. Soc. 41(1), 701 [Q18.8] (1996). 
[12] M. C. Cross and Y. Tu, Phys. Rev. Lett. 75, 834 (1995). 

[13] J. Swift and P. C. Hohenberg, Phys. Rev. A 15, 319 (1977); E. D. Siggia and A. Zip- 
pelius, Phys. Rev. Lett. 47, 835 (1981); P. Mannevill, J. Phys. (Paris) 44, 759 (1983). 

[14] M. C. Cross and A. C. Newell, Physica D 10, 299 (19984); 

[15] Y. Pomeau and P. Manneville, J. Phys. (Paris) 40, L609 (1979); ibid 42, 1067 (1981). 

[16] Y. Pomeau, S. Zaleski and P. Manneville, Phys. Rev. A 27, 2710 (1983); G. Tesauro 
and M. C. Cross, ibid 34, 1363 (1986). 

10 



[17] H. S. Greenside and M. C. Cross, Phys. Rev. A 31, 2492 (1985). 
[18] S. Koga, Prog. Theor. Phys. 67, 164 (1982). 

[19] Following Pomeau and Manneville JOSH , an alternative method is to impose a small 
fluctuation <f)(x,y,t) to the phase 9 m in Eqs. ([I]) and (^), and match the gradients V n 
order by order. This leads to the same results as cited in the text. 

[20] It is straight-forward to extend our calculation to the non-Oberbeck-Boussinesq case in 
which a —gi^ 2 term will be added to Eq. ([!]). Consequently, in addition to mode ([|), 
two extra modes |(An,o + A mt 2e 2t9m + c.c.) are needed. Since these two extra modes are 
generated by the if) 2 term, their amplitudes are both of order A 2 m . 



11 



FIGURES 




FIG. 1. Three-armed global spiral with dislocations at R3, R2 and R\ (R3 < R2 < R\). The 
three arms are represented by the solid, dashed and dotted lines. 



12 



TABLES 

TABLE I. Measured values of uj m R m of global spirals with m = 1, 2, 3 from numerical solu- 
tions. Global spirals are set up with three different sizes R m /L = 0.45, 0.55, 0.65 where L = 32n 
is the radius of the cell. The wave number q m is calculated via Eq. (9), which is within the 
uncertainty of direct measurements q m = 1.00 ± 0.05. 



m 


i? m /32^ 




q m 


1 


0.45 


0.0750 


0.981 




0.55 


0.0742 


0.981 




0.65 


0.0734 


0.981 


2 


0.45 


0.0619 


0.967 




0.55 


0.0630 


0.967 




0.65 


0.0658 


0.965 


3 


0.45 


0.0574 


0.954 




0.55 


0.0568 


0.954 




0.65 


0.0556 


0.955 
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